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ABSTRACT 

We have obtained limited photometric and spectroscopic data for Supernova (SN) 2008bk in NGC 
7793, primarily at > 150 d after explosion. We find that it is a Type Il-Plateau (II-P) SN which most 
closely resembles the low-luminosity SN 1999br in NGC 4900. Given the overall similarity between the 
observed light curves and colors of SNe 2008bk and 1999br, we infer that the total visual extinction 
to SN 2008bk (Ay = 0.065 mag) must be almost entirely due to the Galactic foreground, similar to 
what has been assumed for SN 1999br. We confirm the identification of the putative red supergiant 
(RSG) progenitor star of the SN in high-quality g'r'i' images we had obtained in 2007 at the Gemini- 
South 8 m telescope. Little ambiguity exists in this progenitor identification, qualifying it as the best 
example to date, next to the identification of the star Sk —69° 202 as the progenitor of SN 1987A. 
From a combination of photometry of the Gemini images with that of archival, pre-SN, Very Large 
Telescope JHK S images, we derive an accurate observed spectral energy distribution (SED) for the 
progenitor. We find from nebular strong-intensity emission-line indices for several Hn regions near 
the SN that the metallicity in the environment is likely subsolar (Z ss 0.6Z Q ). The observed SED 
of the star agrees quite well with synthetic SEDs obtained from model red supergiant atmospheres 
with effective temperature T e fj = 3600 ± 50 K. We find, therefore, that the star had a bolometric 
luminosity with respect to the Sun of logtXboi/L©) = 4.57 ± 0.06 and radius i?* = 496 ± 34 R© at 
~ 6 months prior to explosion. Comparing the progenitor's properties with theoretical massive-star 
evolutionary models, we conclude that the RSG progenitor had an initial mass in the range of 8-8.5 
M©. This mass is consistent with, albeit at the low end of, the inferred range of initial masses for 
SN II-P progenitors. It is also consistent with the estimated upper limit on the initial mass of the 
progenitor of SN 1999br, and it agrees with the low initial masses found for the RSG progenitors of 
other low-luminosity SNe II-P. 

Subject headings: supernovae: general — supernovae: individual (SN 2008bk) — stars: late-type — 
stars: evolution — stars: fundamental parameters: other — galaxies: individual 
(NGC 7793) 



1. INTRODUCTION 

The evolution of massive stars, particularly their end- 
points, remains a fascinating, and still unresolved, set of 
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questions and problems that have ramifications in var- 
ious aspects of astrophysics, such as galactic evolution 
and chemical enrichment. Theoretically, we expect stars 
with initial masses greater than 8-9 Mq to terminate cen- 
tral nuclear burning while in the red supergiant (RSG) 
phase of their evolution, un dergo core collapse, and 
explode as super novae (SNe; iWooslev fc Weaver] 119861 : 
IHeger et al.l I2003D . From observations of the progeni- 
tor stars of Type Il-Plateau (II-P) SNe, we find that 
this theoretic al lower mass l imit i s essentially what is be- 
ing inferred (jSmartt et al.l I2009D . The SNe II-P from 
detected stars near this limit have tended to be of 
relatively low luminosity, such as SN 2005cs in M51 
(jMaund. Smartt. fc Danzigerl [2001 ILi et all I2006D and 
SN 2009md in NGC 3389 (jFraser et all 1201 ID . This 
stands in contrast to models which argue that low- 
luminosity SNe II-P can be explained by fallback of ejecta 
on a newly formed black hole, as a result of the core col- 
lapse of a highly massive (> 20 M^) p r ogenitor (e.g., 
i Turatto et alTll998b IZampieri et al.l 120031 : iNomoto et alJ 
I2006I K 

In recent years, SN progenitors have been directly 
identified in pre-SN archival images obtained by the 
Hubble Space Telescope (HST). Of course, direct 
progenitor identification has also been made from the 
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ground prior to that, the most famous example being 
the discovery of the star Sk — 69°202 as the progenitor 
of SN 1987A in t he La rg e Magellanic Cloud (LMC; e.g. 
Gilmozzi et al.l Il987t iSonneborn. Altner. fc Kirshnerl 
19871) . Other cases include t he progenitors of SN 
1961V in NGC 1058 (IBertolal fl96l IZwickvl H961 . 
SN 1978K in NGC 1313 (iRvder et al.l I1993D. and SN 
1993 J in M81 (Aldcring, Humphreys, & Richmond 



Tiil iCohen, Darling, fc Porterl 119951: iVan Dyk et al 
2002( 1 . However, all of these examples are SNe which 
in some way have been peculiar or otherwise unusual. 
The nature of SN 1961V a s a true SN continues to be 
debated (see I Van Dvkll2005l and references therein: also , 
Smith et all 120111: iKochanek. Szczygiel. fc Stanek! 120111: 
Van Dyk fc Mathesonll2011l) . 

SN 2008bk in the nearby, nearly face-on spiral galaxy 
NGC 7793 (see Figured]) has also had its progenitor star 
ide ntified in groun d-based data. The SN was discovered 
bv iMonardl (|2008D on March 25 . 14 (U T dates are used 
throughout this paper). iLi et al.l (|2008f ) were the first to 
attempt to identify the progenitor star as a likely RSG in 
deep archival ground-based BVI images obtained in 2001 
with one of the 8.2 m Unit Telescopes (UT) of the Euro- 
pean Southern Observatory (ESO) Very Large Telescope 
(VLT). They used the Monard discovery images and also 
images taken on March 28.41 with the 0.41 m PROMPT 
2 telescope at the Cerro Tololo Inter- American Observa- 
tory (CTIO) as references. This search for the progenitor 
was attempted even before the SN had been classified. 
However, Li et al. concluded that SN 2008bk was likely 
a core-collapse SN, since a star was detected at the loca- 
tion; it is far less likely that we would be able to detect 
the progenitor, or the companion to the progenitor, of a 
thermonu clear Type la SN outside of the Local Group 
(see, e.g.. lLi et al.ll2011l). The S N was subsequently clas- 
sified by iMorrell fc Stritzingerl 1)20081 ) as Type II with 
an age of 36 d after explosion on April 12.4, based on 
a comparison with the well-studied SN II-P 1999em in 
N GC 1637. 

iMaoz fc Mannuccl (120081 ) detected and measured mag- 
nitudes for the progenitor star in archival near-infrared 
(IR) J- and JC, s -band images, also obtained with the VLT 
in 2005. IMattila et all (120081 ) further analyzed the VLT 
optical images, as well as VLT ii-band data obtained in 
2007, and combined these pre-SN results with precise as- 
trometry provided by post-SN K s VLT adaptive optics 
imaging. From their analysis they concluded that the 
total extinction toward the progenitor is Ay — 1.0 ± 0.5 
mag, and that the star therefore was likely of type M4 I, 
with an initial mass of 8.5 ±1.0 M©. 

In this paper, with limited photometry and spec- 
troscopy of the SN itself, we show that SN 2008bk is 
a SN II-P very similar to the low-lum inosity SN 1999br 
in NGC 4900 (jPastorello et al.l 120041 ). Furthermore, we 
measure accurate photometry for the RSG progenitor 
from superior-quality g'r'i' images obtained in 2007 with 
the Gemini-South 8 m telescope, as well as from the 
archival, pre-SN, near-IR VLT images. We compare 
this photometry to spectral energy distributions (SEDs) 
obtained from stellar atmosphere models for RSGs, to 
infer the star's effective temperature. Together with 
the inferred luminosity for the star, we compare these 
properties to recent theoretical models for massive-star 
evolution, to estimate the initial mass for the progeni- 



tor. We also compa re our results to those obtained by 
IMattila et all (|2008l ). 

2. SUPERNOVA OBSERVATIONS 

Although highly incomplete in terms of photomet- 
ric and spectroscopic coverage, our observations of SN 
2008bk are still sufficient to determine its general prop- 
erties and to allow for comparison with other SNe II. A 
more complete analysis of the properties of SN 2008bk is 
forthcoming (Pignata et al., in prep.). 

2.1. Photometry 

We obtained BVRI images of SN 2008bk, at eight 
epochs late in its evolution: four at the Calar Alto 2.2 
m telescope + the CAFOS SiTe (scale 0'/53 pix~\ field 
of view 9' x 9'), two at the 2.5 m Nordic Optical Tele- 
scope (NOT; Roque de los Muchachos Observatory, La 
Palma, Spain) + ALFOSC (0'.'19 pix" 1 , field of view 
6f5 x 6.'5), and two using the SMARTS 1.3 m telescope 
+ ANDICAM at CTIO. These images were processed 
using standard IRAfF^l procedures for image trimming, 
bias subtraction, and flat-field corrections. Additionally, 
we add several relatively early-time epochs at V , with an 
image obtained at the Las Campanas Observatory (LCO) 
2.5 m duPont telescope on 2008 April 13 and six images 
obtained with the Global Rent-A-Scope (GRAS) 0.32 m 
and 0.41 m telescopes at epochs spanning the plateau 
to early on the exponential tail. We have also included 
BV data, obtaine d by the CHilean Au tomatic Supernova 
sEarch (CHASE: IPignata et al.ll200l with PROMPT 2. 

We extracted instrumental photometry for the SN and 
reference stars from these images, generally, vi a point- 
sprea d function (PSF) fitting using DAOPHOT (jStetsonl 
1987) within IRAF. When the shapes of the stars in im- 
ages deviated substantially from a "well-behaved" PSF, 
we relied only on aperture photometry. For this rea- 
son, we calibrated the SN and reference star instrumental 
photometry from all of the images to existing photome- 
try for the star fi eld around the host galax y NGC 7793 
was observed by lLarsen fc Richtlerl (119991 ) in UBVRI 
with the Danish 1.54 m telescope + the Danish Faint 
Object Spectrograph and Camera on 1997 September 
6. We obtained these images, posted for public distri- 
bution in FITS format, from the NASA/IPAC Extra- 
galactic Database (NEEF^). We extracted instrumental 
magnitudes from the images using IRAF, first with a 4" 
aperture and then with a PSF for each band. Calibra- 
tion was established using photometry of a number of 
isolated reference stars through a 20" aperture, which 
matches the aperture used to establi sh the calibra tion 
through standard-star observations bv iLarsenl ([19991 his 
Table 1); see Figure [1] and Table Q] for the locations and 
magnitudes of these stars. We also applied the appro- 
priate extinction corrections (S. S. Larsen 2010, private 
communication) . 

We can compare our BVR photometry with those 
magnitudes provided in NOMAEP^I: however, the un- 

12 IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foun- 
dation (NSF). 

13 http:/ /nedwww. ipac.caltech.edu. 

14 http:/ /www. nofs.navy.mil/nomad. html. 
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certainties in these latter magnitudes are substantial, 
resulting in large scatter in the comparison: for B, 
NOMAD-Ours = -0.32 ±0.20; for V, -0.51 ±0.15; and 
for R, — 0.39±0.19. A far better comparison can be made 
at / with the more reliable photometry from the DE- 
NIS database ([DENIS Consortiuml[2005f i: DENIS -Ours 
= 0.02 ± 0.05 mag. The good agreement in the / mag- 
nitudes provides us with some confidence, overall, in the 
calibration of our photometry. 

As many reference stars as possible were generally 
adopted for the calibration; however, for the shallowest 
images, we had to use a combination of the magnitudes 
of star K (the brightest star in Table [1]) and the less 
precise BV photometry from the bright NOMAD stars 
seen in these images. The resulting photometry of SN 
2008bk is given in Table [2j The uncertainties in the SN 
photometry in all four bands arise from the measurement 
uncertainties from DAOPH OT, the uncert ainties in the 
absolute calibration given bv lLarsenl (119991) . and the un- 
certainties in transforming the instrumental magnitudes 
to standard magnitudes, based on the reference stars, all 
added in quadrature. 

The BVRI light curves for SN 2008bk are shown in 
Figure [2] As we have noted, many of the multi-band 
points are beyond the plateau phase and on the expo- 
nential tail. However, the V photometry at earlier times 
provides us with adequate leverage on the overall photo- 
metric evolution of SN 2008bk, from soon after explosion 
through the onset of this tail. The lone early-time B — V 
datapoint provides us with some leverage on the initial 
color of the SN. 

For comparison, we also show in Figure [5] the 
V light c urves for bo t h th e "normal" SN II- 
P 1999em (IHamuv et all 120011 : iLeonard et all l2002t 
lElmhamdi et al l I2003D and the sub luminous SN II-P 
2005cs in M51 (IPastorello et al.ll2009D. and, for all bands , 
the subluminous SN II-P 1999br ([Pastorello et aLll2004D . 
We have attempted to match, by eye, the brightnesses of 
these comparison SNe with SN 2008bk toward the end of 
the plateau phase at V. Although we have very limited 
coverage at the end of the plateau, it appears that this 
phase ended very near JD 2454651 (i.e., 2008 July 3), 
implying that the plateau lasted ~ 100 d after discov- 
ery. (Additionally, the end of the plateau phase for SN 
1999br is not well defined; IPastorello eFaT] 12004 ) The 
SN 2008bk light curves are shown in Figure [2] relative to 
that epoch. 

In Figure [3] we show the absolute V light curve 
for SN 2008bk, which results from correcting for the 
assume d ex tinction (see § 14. ip and distance modulus 
(see § 14. 2|) to the SN. We also show for compari- 
son the absolute V light curves for SN 1999em, cor- 
rected for the tota l extinction, Ay = 0.31 mag (e.g., 
IHamuv et al.l 120011). and the Cepheid distance modulus 
to its host galaxy ([Leonard et a l. 2003); SN 2005 c s, cor - 
rected for the extinction from IPastorello" eFall (2009) 
and the distance modulus from iTonrv et all (f2001): 
and SN 1999br (extin ction and distance modulus from 
IPastorello et al.ll2004D . We estimate that SN 2008bk had 
an absolute magnitude near maximum brightness dur- 
ing the plateau phase of M v = —14.8 mag, which is 
~ 2.1 mag less luminous than SN 1999em. Although 
the plateau brightness for SN 2008bk is very similar to 
the M v = -14.9 mag for the low- luminosity SN II-P 



2005cs (|Pastorello et al.l [20091) . SN 2008bk did not de- 
cline nearly as rapidly, post-maximum, nor is it nearly 
as subluminous on the tail as SN 2005cs. One can see 
from Figure [2] that SN 2008bk exhibits a striking photo- 
metric resemblance to SN 1999br at all bands. However, 
a noticeable offset exists between the absolute bright- 
nesses of SNe 2008bk and 1999br; we comment on this 
more in §[5] We will demonstrate in g2^]that SN 2008bk 
also spectroscopically most closely resembles SN 1999br. 
We therefore consider SN 2008bk to be a low-luminosity 
SN II-P, quite similar to SN 1999br. 

The photometric comparison with SN 1999br (assum- 
ing the explosion da te for SN 1999br is JD 2451278; 
IPastorello et al.ll2004D indicates that the explosion epoch 
for SN 2008bk was ~ JD 2454550, or - 2008 March 24. 
Monard, therefore, likely discovered the SN < 1 d after 
explosion. 

2.2. Spectroscopy 

We obtained several spectra of SN 2008bk at late times, 
using the Blanco 4 m telescope at CTIO with the RC- 
Spectrograph in 2008 July and using the Calar Alto 2.2 
m telescope with the CAFOS SiTe imaging spectrograph 
in 2008 September. The details of these observations 
are given in Table [3] The spectra were reduced using 
standard IRAF routines. The preprocessing of the spec- 
troscopic data (trimming, bias, overscan, and flat-field 
correction) were the same as for the images (see tj2.ip . 
The one-dimensional spectra were wavelength calibrated 
by comparison with arc-lamp spectra obtained during 
the same night and with the same instrumental con- 
figuration, and flux calibrated using spectrophotometric 
standard stars. The zeropoint of the wavelength calibra- 
tion was verified against the bright night-sky emission 
lines. The standard-star spectra were al so used to model 
and r emove the telluric absorption fe.g..lMathcso n~et al.l 
[20001) . 

The resulting spectra are shown in Figure [4] Since the 
spectra obtained at CTIO are only spaced in time each 
by a day, the features in the spectra had varied little 
in that time interval. Thus, we show only a represen- 
tative spectrum from July 24. The same is true for the 
two spectra obtained at Calar Alto in mid-September; 
we display only the spectrum from September 16. The 
July 24 spectrum is ~ 21 d past the plateau phase, while 
the September 16 and 30 spectra are during the nebular 
phase, at ~ 75 and ~ 89 d past the plateau. In Figured] 
we show the SN 2008bk spectrum from Ju ly 24, and also 
spect ra of the normal SN II-P 1999em ([Leonard et al.l 
2002) at 1 d and 29 d past the plateau phase, as well as 
a spe ctrum of the sublumin ous SN II-P 1999br in NGC 
4900 ([Pastorello et alJl200l ) at essentially the end of the 
plateau. The spectra for SN 1999em were obtained from 
the SUSPECT SN database^. A. Pastorello kindly pro- 
vided the spectrum from the Padova SN data archive of 
SN 1999br from 1999 July 20. 

Both the strength and width of the various spectral 
lines, and the overall continuum shape, for SN 2008bk 
compare rather poorly with SN 1999em. A far bet- 
ter comparison, overall, can be made with SN 1999br. 
The relative strength of, for example, the [Call] AA7291, 
7324 feature in the SN 2008bk spectrum, compared with 

5 http:/ /suspect. nhn.ou.edu/~suspect. 
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that in the SN 1999br spectrum, again indicate that SN 
2008bk is also a low-luminosity SN II-P. 

3. DETECTION OF THE PROGENITOR STAR 

ILi et al.l (]2008[ ) measured a precise position for the SN 
of a(J2000) = 23 h 57 m 50M2, <5(J2000) = -32°33'21'/5. 
Within the uncertainty of this position (±0'.'2) they iso- 
lated a star in the VLT /-band image at this location. 
The star is not detected in the B and V images from 
the same epoch. This candidate progenitor lies within 
a prominent cluster of bright blue and red supergiant 
stars. From its environment and its lack of detection in 
the bluer bands, Li et al. conjectured that the progenitor 
was a RSG. 

We have since improved the astrometric agreement, 
with an uncertainty of 0'.'13, between the SN and the star 
with a V image of the SN from the LCO 2.5 m duPont 
telescope, obtained on 2008 April 13. We have refined 
the absolute position for the SN to end figures 50?50 and 
20"7 (with uncertainty ±0'.'2). 

We obtained deep g'r'i' images in 2007 September and 
October with GMOS (Crampton et al. 2000) on Gemini- 
South, to study the outer disk of NGC 7793. The im- 
ages were obtained in queue observing mode as part of 
program GS2007B-Q-57 (PI: Davidge). The detector in 
GMOS is a mosaic of three CCDs with a raw image scale 
of 0'.'0727 pix" 1 . The data were binned 2x2 pixels 2 
to better match the delivered image quality. Eight 200- 
s exposures were obtained in g\ and eight 300-s expo- 
sures were obtained in both r' and i 1 . The raw images 
in each band were bias subtracted, divided by a twilight 
flat, and then stacked. The progenitor star and its envi- 
ronment are shown in Figure [6j From these images we 
measure an absolute position for the star to end figures 
50! 48 and 20'/9 (±072). To within the absolute astromet- 
ric uncertainties, the position of the identified progenitor 
star agrees very well with the SN position (see above). 

We extracted instrumental photometry from the g'r'i' 
images via PSF fitting within IRAF/DAOPHOT. Since 
the Gemini program was intended to be carried out dur- 
ing nonphotometric conditions as part of the so-called 
"Poor Weather Queue," this photometry must be cal- 
ibrated using other observations. Following Welch et 
al. (2007), we then transformed the g'r'i' instrumental 
magnitudes from the Gemini images to the calibrated 
VRI magnitudes from the Danish telescope images (see 
§2.1j) . using a number of bright (V — 16.1-20.4 mag), iso- 
lated stars common to both sets of images. The resulting 
transformations are 

V = 0.7803.g' - 0.0181r' + 0.2768i' + zeropoint, 

R = O.OOOOg' + 0.7762r' + 0.2419?' + zeropoint, 

I = 0.0000.g' - 0.0722r' + 1.1055i' + zeropoint. 

We note that the coefficients in the transformation above 
are similar to the ones derived by Welch et al. as ap- 
plied to a different scientific context. The uncertainties 
in these transformations are 0.08 mag for V, 0.05 mag 
for R, and 0.07 mag for /. 

The resulting VRI magnitudes for the SN 2008bk pro- 
genitor from the Gemini images are given in Table H] The 
uncertainties in the magnitudes include the measurement 
uncertainties from DAOPHOT, the uncertainties in the 



calibration of the Danish telescope images of the galaxy 
(|Larsen|[l999l his Table 1), and the uncertainties in the 
transformations above, all added in quadrature. 

We also reanalyzed and remeasured the photometry 
from the near-IR VLT IS AAC and HAWK-I JHK S im- 
ages, initially analyz ed bv lMaoz fc Mannu cci (2008]) and 
iMattila et al.l (|2008f) . Photometry was performed, again 
using PSF fitting with IRAF/DAOPHOT, and the in- 
strumental magnitudes were calibrated to 2MASS stars 
in the fields. The resulting photometry is given in Ta- 
ble |U The uncertainties in the magnitudes include the 
measurement uncertainties from DAOPHOT, the uncer- 
tainties in the 2MASS star magnitudes, and the uncer- 
tainties in the transformations using 2MASS. Later (see 
§ 14.41). we will need these near-IR magnitudes in the 
iBessell fc Brett] (|1988l ) system, particularly for K. For 
this reaso n, we use the tran sformations from the 2MASS 
system in iCarpenterl (|2001b his Appendix A) to obtain 
j = 19.34, H = 18.56, and K = 18.18 mag (formal un- 
certainties in these transformations are all < 0.01 mag). 

4. PROPERTIES OF THE PROGENITOR STAR AND ITS 
ENVIRONMENT 

We will now place the SN 2008bk progenitor in a 
Hertzsprung- Russell diagram (HRD), to attempt to esti- 
mate its initial mass. Whereas the photometry of the 
progenitor star has relatively small uncertainties, the 
next steps, discussed below, involved in analyzing this 
photometry possess significantly larger uncertainties. 

4.1. Extinction to the Supernova 

Here we estimate the extinction to the SN. The 
startlingly good agreement, seen in Figure [2] between 
the light curves of both SN 2008bk and SN 1999br 
in all observed bands (the SN 1999br curves have all 
been uniformly adjusted by —4.6 mag to match, by 
eye, those of SN 2008bk) implies that the extinction to- 
ward SN 2008bk is quite similar to that of SN 1999br. 
iSmartt et al.l (|2009f) assume only the Galactic foreground 
contribution to the extinction toward SN 1999br: Ay = 
0.065 mag. We can then logically infer that the extinc- 
tion toward SN 2008bk is comparably quite low; we as- 
sume Ay = 0.065 mag. Interestingly, the estimate for 
the Galactic line-of-sight extinction toward the SN is, 
in fac t, Ay — 0.065 mag (jSchlegel. Finkbeiner. fc Davisl 
Il998h . which would imply that the extinction in the im- 
mediate galactic environment of SN 2008bk is negligible. 
From a preliminary analysis of the photospheric spectra 
for SN 2008bk, no sign exists of Na I D absorption due 
to the host galaxy. Additionally, the average internal ex- 
tincti on in NGC 7793 is g enerally quite low (Ab = 0.12 
mag: iPierce fc Tullvlll992 ). The low extinction value is 
consistent with the lack of emission at 8.0 /jm (normally 
attributed to emission from polycyclic aromatic hydro- 
carbons in interstellar dust clouds) detected at the posi- 
tion of SN 2008bk in archival images of NGC 7793 ob- 
tained using the InfraRed Array Camera (IRAC) on the 
Spitzer Space Telescope. It is also consistent with the 
fact that the SN is within what appears to be a nebular, 
possibly evacuated, bubble (see Figure [7]). 

4.2. Distance Modulus to the Host Galaxy 

We also requ ire an estimate of t he di stance modu- 
lus to the SN. iPuche fc Carignanl (|1988t ) calculated a 
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mean distance modulus, 27.64 ± 0.19 mag, to NGC 7793 
from several dis tance indicators, includin g the Tully- 
Fisher relation. Kara chentsev et all (|2003| ) measured a 
distance modulus of 27.96 ± 0.24 Mpc (somewhat more 
distant, but in good agreement with the Puche & Carig- 
nan value) , using HST images of the host galaxy to esti- 
mate the "tip of the red giant branch" (TRGB) bright- 
ness to be I = 23.95 ± 0.22 mag. (These HST images, 
unfortunatel y, did not cover the SN site; furthermore, 
we note that iMould fc SakaH^OOa using the TRG B cal- 
ibration of Mf = -4.05 mag bv lRizzi et al.l [20071 found 
/i° = 27.78 mag; however, their assumed /-band extinc- 
tion, Aj = 0.22 mag, is almost certainly too large.) 

The distance modulus has been re cently determined 
using Cepheids in the host galaxy by iPietrzvhski et all 
(20101) to be /I = 27.68 ± 0.05 mag (random) ±0.08 
mag (systematic). Those authors note that their dis- 
tance modulus also agrees very well with a recent TRGB 
measu rement, ^° = 27.72 ± 0.08 mag, by Uacobs et alJ 
(2009). Hereafter, we adopt the Cepheid distance mod- 
ulus and, conservatively, assume the (larger) systematic 
uncertainty in that estimate. 

4.3. Metallicity of the Supernova Environment 

Furthermore, we need to know the metallicity in the SN 
environment, since we are comparing the progenitor's ob- 
served characteristics with estimates from models which 
are generated assuming a particular metallicity. 

We consider the metallicity to be equivalent to the local 
oxygen abundance. The most accurate estimate of the 
metallicity ideally would be from the immediate SN en- 
vironment itself, or failing that, from regions of the host 
galaxy in the general vicinity of the SN. We show sev- 
eral such H II regions in Figure which is a continuum- 
subtracted Ha (+ [Nil]) image of NGC 7793, released as 
an ancillar y data product by the Spitzer SINGS Legacy 
project (jKennicutt et al.l [20031) . In the spectroscopic 
observations of the SN obtained on 2008 September 30 
at Calar Alto, light from two Hn regions, "A" and "B," 
approximately to the south of the SN also fell within the 
C AFOS slit. 

iPettini fe Page! (|2004D calibrate relations between 
the O abundance and the strong-inten sity line indices 
N 2 = l og [Nil] A6584/Ha (see, e.g.. iTerlevich et all 
199U IStorchi-Bergmann, Calzetti. fc Kinney! 119941: 
Denicolo. Terlevich. fc Terlevichll2002t) and 03N2 = log 
(([Oiii] A5007/H/3)/([Nn] A6584/Ha)). These two 
indices are relatively unaffected by extinction, since 
they depend on ratios of intensities for two spectral 
lines essentially contiguous in wavelength. The la 
uncertainties in the O abundance from t hese two indices 
are ± 0.18 and ±0.14 dex, respectively ([Pettini fc Pagell 
2004). We measured the intensities for these various 
emission lines, using standard routines in IRAF, and 
estimate these line indices for these two H II regions. 

The resulting values of the indices and the computed 
abundances are given in Table O We estimate that 
our measurement errors lead to uncertainties in the 
abundances of, conservatively, ±0.05 dex. Addition- 
ally, we can co nsider the measurements for the H II re- 
gions #27 fromlBibbv fc Crqwtherl (|2010h and W13 from 
iMcCall. Rvbski. fc Shields! (|1985[ ). These two regions 

16 Available at http://irsa.ipac.caltech.edu 



are labeled in Figure [7] and listed in Table [SJ We also 
include in the table the deprojected angular distance 
of each Hn region from the SN, assuming the values 
of the host-galaxy inclina tion and position angle from 
IBibbv fc Crowtherl (|2010D and references therein. The 
linear distances in the table are assuming a distance to 
the host galaxy of 3.4 Mpc (see 

One can see that the O abundances of these H II re- 
gions are all consistent with an approximate range of 8.4 
to 8.5 dex. The region nearest to (~ 200 pc from) the 
SN, in particular, has 12 + log(0/H) = 8.45. Assum- 
ing the solar p hotospheric abundanc e is 12 + log(0/H) 
= 8.66 ± 0.05 ([Asplund et al.ll2004D . we can infer that 
the metallicity in the SN environment is « 0.6 solar. 
The mean O abundance of H II regions in the LMC i s 12 
± log(0/H) = 8.37 ±0.22 ([Russell fc Dopital [1991 . or 
~ 0.5 solar. Therefore, the metallicity at the SN site is 
most likely not solar, but also is likely not quite as low 
as the mean LMC nebular met allicity. Howeve r , with 
the level of uncertainties in the IPettini fc Pagell (2004) 
calibration relations, in our measurements, and in the 
assumed solar and LMC O abundances, we analyze the 
progenitor photometry considering a range of metallici- 
ties. We will show that the inferred intrinsic properties 
for the progenitor, nonetheless, are most consistent with 
subsolar metallicity in the SN environment. 

4.4. Initial Mass of the Progenitor 

We have assumed a relatively low extinction toward 
SN 2008bk, based on the comparison with SN 1999br 
f q2.1[) . We further assume that this extinction f t|4.1l) is 
also appropriate to apply to the progenitor star. From 
the Gemini optical data, this star, therefore, had an ab- 
solute magnitude M v = —4.94 and a dereddened color 
(V — I)o « 2.19 mag, which imply that the progeni- 
tor w as a lower- luminosity (class lb ), early M-type RSG 
(e.g., Humphreys fc McElrovlll984D . 

From the observed optical and near-IR photometry 
in Table |U we show the SED of the star in Figure [3 
For comparison we show synthetic SEDs, generated us- 
ing SYNPHOT ap plied to the MAR CSEZl model stellar 
atmospheres {Custafss on et al.l [2008) for RSGs with ef- 
fective temperature T c g = 3600 K and surface gravity 
log<? = 0.0, assuming spherical geometry and subsolar 
metallicity ([M/H] = -0.25, where [M/H] denotes the 
ratio of the general metal abundance to H, relative to 
the Sun). We also show the models at solar metallicity 
([M/H] = 0) for this temperature, and one can see little 
difference with the subsolar models. The MARCS mod- 
els are for stars with initial masses 5 and 15 M Q (we will 
show that the initial mass of the SN 2008bk progenitor 
was in between these two values) and have microturbu- 
lence velocities spanning from 1 to 5 km s~ 1 . The models 
were all reddened b y the assumed value for the p rogen - 
itor, following the ICardelli. Clayton, fc Mathisl (|1989f) 
reddening law, and were normalized to the observed V 
magnitude. The agreement of all of these models with 
the observed SED is quite good, to within the uncertain- 
ties in the photometry of the progenitor. (However, the 
models do diverge by > la from the observation at R, 
and compare less well with the observations at H.) We 
therefore consider the effective temperature of the star 

17 http://marcs.astro.uu.se/. 
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to be T e g = 3600 K, with a conservative uncertainty of 
±50 K. 

This effective temperature c orresponds to spectra l type 
M2-M2.5 at LMC metallicity ([Levesque et al.ll2006t their 
Table 4), and M3 at solar metallicity ([Levesque et abl 
[2001 their Table 5). The metallicity in the SN 2008bk 
environment is li kely s omewhere in between f§ l4.3p . From 
iLevesaue et all ([20051 ; their Table 6) for Galactic (solar 
metallicity) RSGs, the V and K bolometric corrections 
at T cS = 3600 K are BCy = -1.75 and BG K = 2.84 
mag, respectively. At LMC metallicity, the bolometric 
corrections vary little from the solar val ues at this T e g: 
BCy = -1.78 and BC K = 2.84 mag (jLevesque et al.l 
2006). The absolute bolometric magnitude of the star, 
therefore, from V is Mboi = —6.70 ±0.12 mag, and from 
K it is Mboi = —6.67 ±0.13 mag. The total uncertainties 
arise from the uncertainties in the measured photometry 
for the star and in the distance modulus (see H4.2p . and 
at V, in the range in BCy- That these two estimates 
of Mboi, established from the two different photometric 
measurements, agree so well to within the uncertainties, 
instills us with considerable confidence in the star's pho- 
tometry and in our estimates of the extinction and the 
star's effective temperature. Assuming Mboi(©) = 4.74 
mag, the bolometric luminosity relative to the Sun is then 
log(L b oi/L Q ) = 4.57 ±0.06 from both V and K. The lu- 
minosity and effective temperature together indicate that 
the RSG's radius, at ~ 6 months prior to explosion, is 
R* = 496 ±34 R Q . 

The inferred intrinsic properties of the progenitor of 
SN 2008bk are shown in the HRD in Figure M We can 
also compare the star's locus in the HRD with massive- 
star evolutionary tracks of various initial masses, and, 
given the uncertainties mentioned in § I4.3| we take into 
consideration tracks at several metallicities. In panel 
(a) of the figure, we compare to tracks at solar metal- 
licity _XZ_^_0J)2VU„the Cambridge STARS tracks 
fromlEldridge fc Toutl (120041) and t he Geneva tracks from 
iHirschi. Mevnet. fc Maederl ([2004D . which include either 
stellar equatorial rotation, at v- m i — 300 km s" 1 , or no 
rotation. One can see that these tracks would imply that 
the initial mass of the progenitor is in the range of ~ 10- 
12 Mq. However, these tracks all terminate at cooler ef- 
fective temperatures and higher luminosities (and larger 
stellar radii) than is the case for the progenitor. In total, 
the solar metallicity tracks compare rather poorly with 
the progenitor's inferred properties, and we consider it 
unlikely that the larger initial masses implied by these 
tracks apply to the SN 2008bk progenitor. 

In Figure [9j panel (b), we show a comparison with 
the STARS tracks at Z = 0.01 (i.e., at Z/Z© = 0.5; 
the Geneva tracks have not been published at subso- 
lar metallicity in this mass range). The 8 and 8.5 M© 
model tracks both terminate with T e g = 3567 K, and 
with log(Lboi/L©) = 4.59 and 4.63 (respectively) and 
R± = 514 and 541 R© (respectively), which are essen- 
tially consistent, to within the uncertainties, with those 
properties of the progenitor. Note that the 9 M© model 
terminates at a luminosity that is higher than the range 
we have inferred for the progenitor's luminosity. Note 
also that the 7.5 M© track terminates at a much cooler 
temperature and much higher luminosity, and therefore 
we can discount it entirely. In panel (c) , we compare with 



the STARS models for Z = 0.008 (i.e., Z/Z© = 0.4), al- 
though it is less likely that the metallicity in the SN 
environment is quite this low (§ I4.3|) . The 8 and 8.5 
M© tracks terminate with T c g = 3615 and 3622 K, 
log(L b oi/L©) = 4.60 and 4.63, and R* = 508 and 527 R© 
(respectively), which again agree, to within the uncer- 
tainties, with the star's locus in the HRD (the exception 
being the 7.5 M© track, which terminates at a stage too 
cool and too luminous than what is observed, and again, 
we can discount it). Much better agreement exists be- 
tween the star's properties and the subsolar metallicity 
tracks, which, all told, imply that the star's initial mass 
was ~ 8-8.5 M©. These stellar evolutionary models all 
have lost ^0.3 M© of mass prior to reaching their ter- 
mini. 

4.5. Comparison with Mat tila et al\ I 20 OA ) 

We can make a direct comparison of o ur results for 
the SN 2008bk progenitor with those of IMattila et al.l 
(2008). One can see from Figure [10] that all of our 
magnitude measurements for the star a r e > l a brighter 
than those measured by IMattila et al.1 (|2008f ) , particu- 
larly their measurement at / and their upper limit at 
V (which is substantially fainter than our actual detec- 
tion in this band). We cannot provide an explanation 
for this difference. However, we feel confident in the 
veracity and self-consistency of our photometry, based 
on comparison with pre-existing photometry of stars in 
the SN field (§ I2.1j) . and on the relatively low disper- 
sions in the transformation of g'r'i' photometry to VRI 
from the Gemini images and of the near-IR instrumen- 
tal photometry to 2MASS magnitudes from the VLT im- 
ages ( §[3l). Additionally, the Elias. Frogel. fc Humphreys! 
(1985) RSG intrinsic colors, specifically (V—K) , are sys- 
tematically blu er, for a given T P a and spectral type, than 
those found bv lLevesque et al. (2006). The net effect, as 
seen in Figure fTUl is that higher extinction and reddening 
are required to allow for a comparison (and not a particu- 
larly good one!) of the model SEP for a M4I supergiant 
from [Elias. Fr ogel. fc Humphrey s (1985) with the pho- 
tometry from IMattila et all ([20081 ). Simply for the sake 
of argument, we also show in the figure the SED for a 
RSG from the MARCS stellar atmos phere models with 
T c ff = 3500 K (which Mattila et al.l assume to be the 
equivalent of a M4I star at LMC metallicity), adopting 
our extinction to the SN. Th is model actua l ly pro vides 
a better comparison with the IMattila et al.l (120081) pho - 
tometry than do the lElias. Frogel. fc Humph reysl ([19851 ) 
intrinsic colors, and this not only would imply that the 
effective temperature is cooler than what we have found, 
but does so in agree ment with a fa r lowe r value for the 
extinction than what IMattila et al.l (|2008[ ) had assumed. 
The implied cooler temperature is driven primarily by 
the significantly fainter /-band magnitude and uppe r 
limit to the V brightness found by IMattila et al.l ([2 0081. 

The fainter K-band measurement by IMattila et al.l 
(2008), combined with the implied cooler effective tem- 
perature, which corre s ponds to a larger (positive) BCk 
([Levesque et al.l [200 5: 2006), as well as the lower value 
for extinction, all lead to a much lower luminosity for the 
SN progenitor than what we have found to be the case. A 
comparison of the luminosity and effective temperature 
with massive-star evolutionary models in a HRD shown 
in Figure [TT] implies that the SN progenitor would have 
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an initial mass well below the threshold mass necessary 
for core collapse. It is only by virtue of their assumption 
that Ay ~ 1.0 mag, which, based on their photometry, 
results in the star being mo re luminous than what we 
find here, that iMattila et al.l arrive at its initial mass be- 
ing Mi n i ~ 8.5 M Q . As we presented in § 14. 1| we have 
found no evidence for such a large extinction to the SN. 

Estimating the reddening to the progenitor through 
comparison with assumed intrinsic properties for 
RSGs is almost ce rtainly fraught with uncertainties. 
IMattila et al.1 (|2008l ) also attempted to estimate the ex- 
tinction bas ed on the Balmer decrement for an H II re- 
gion (W13; IMcCall. Rvbski. fc Shields! H985h . which is 
1.'5 from the SN 2008bk site. In general, it is far more 
accurate to estimate the reddening from the SN light 
curves and colors, particularly during the plateau phase, 
where the colors of SNe II- P are remarkably similar (e.g., 
IBersten fc Hamuvl 120091 ). Although our light curves for 
SN 2008bk are very limited in coverage, we therefore be- 
lieve that the comparison with the more extensive SN 
1999br curves provide a reliable estimate of the extinc- 
tion to SN 2008bk. 

We note that IMattila et all (12008ft assume d the dis- 
tance modulus from lKarachentsev et al.l (|2003l ); however, 
this contributes relatively little to the overall difference 
between our two studies. Furthermore, our assessments 
of the meta llicity in the S N 20 08bk environment dif- 
fer as well. IMattila et al.l (120081) estimated the metal- 
licity, assuming the IPilvugin. Vflchez. fc Continil ([200l 
metallicity-radial distance relation for NGC 7793, and 
concluded tha t the O abundance was 12 + log(0/H) 
= 8.2 ± 0.1. iBibbv fc Crowtherl ([20101) also have de- 
termined the metallicity gradient in NGC 7793. If 
we assume the deprojected radial d istance, r, of SN 
2008bk given by IMattila et al.l (|2008l 3'A7) and the ra- 
dial distance at B = 2 5 mag arcsec -2 , R25 — 4f65 
(jBibbv fc Crowtherl 120101 and references therein), then 
for r/i?25 = 0.75, the abundance is a slightly higher 12 
+ log(0/H) = 8.34 ±0.05. We have shown in § S3] from 
direct estimation of the O abundance in neighboring H 11 
regions that 12 + log(0/H) w 8.4-8.5. 

5. DISCUSSION AND CONCLUSIONS 

We conclude that SN 2008bk in NGC 7793 is most 
likely a low-luminosity SN II-P, similar in characteristics 
to SN 1999br. We have also confirmed the identifica- 
tion of an RSG as the SN progenitor. This star is read- 
ily detected in prc-SN, high-quality optical images from 
the Gemini-South Observatory. Although SN 1987A is 
closer, its association with the star Sk — 69°2 02 certain 
its m ulti-band brightness and color known (Isserstedt 
119751) . and, more i mportantly, its spect ral type un- 
ambiguously defined ([Rousseau et allll978D . SN 2008bk 
ranks as the "second best-characterized" SN progenitor 
identification to date of any SN type. For this reason, its 
observed properties are well constrained. 

We have assumed, based on the photometric compar- 
ison with SN 1999br, that the extinction toward SN 
2008bk is quite low, Ay — 0.065 mag, entirely consis- 
tent with Galactic foreground reddening. We have also 
estimated the metallicity in the SN environment from 
neighboring Hn regions and find that it is consistent 
with subsolar metallicity (Z w 0.6Z Q ). Furthermore, 
based on comparison of the observed SED with recent 



models for RSG stellar atmospheres, we are able to as- 
sign an effective temperature of T e g = 3600 ± 50 K, and 
an absolute bolometric luminosity (relative to the Sun) of 
log(£boi/L©) = 4.57 ± 0.06, both quantities with rather 
small uncertainties. We note the very good agreement in 
the bolometric magnitudes derived independently from 
both the observed V and K. Through comparison with 
recent theoretical massive-star evolutionary models for 
subsolar metallicity, we conclude that the progenitor of 
SN 2008bk had an initial mass in the range of 8-8.5 M Q . 
We note that this inferred mass range is consistent with 
the upper limit derived for th e initial mass of the progen- 
itor of SN 1999br (< 15 M ; ISmartt et~aTl [20091 ) . and is 
also within the initial mass range for RSG progenitors of 
SNe II-P, M ini = 8.5±^ 5 to 16.5 ± 1.5 M Q (ISmartt et all 
[2001 . 

Further constraints on the initial mass of the 
SN 2008bk progenitor can be obtained through 
comparison with recent theoretical, radiative- 
transfer/hydrodynamical simulations of SNe II-P. 
We can compare our estimate of the star's radius 
near explosion, i?* = 496 ± 34 R , to those input to 
the pre-SN models generated at s olar metallicity by 
iDessart. Livne. fc Waldmanl ((2010b): their nonrotating 
11 M Q model has i?* = 587 R , while their rotating 
10 M has R* = 552 R . This implies that the SN 
2008bk progenitor had to be < 10 M , consistent 
with what we infer from comparison with theoretical 
massive-star evolutionary tracks. 

The low tail-phase luminosity for SNe II of the SN 
1999br variety has been explained by a low 56 Ni mass 
produced in t he ex plosion of a lower mass star. As 
ISmartt eF al. (2009|) point out, nuclear statistical equi- 
librium is therefore reached only in a thin O-Si-rich layer 
around the core in such stars. Low progenitor masses are 
consistent with what has been directly found for other 
low-lumi nosity SNe II-P, such as SN 1999br, SN 2003gd 
in M74 dVan Dvk. Li. fc Filippenkol l2003USmartt et all 
200llHendrv et al.ll2005D . and SN 2005cs (|Li et alll200a 
Pastorello et al.l 120091 ). The lower initial mass estimate 
for the progenitor of SN 2008bk presented here is con- 
sistent with this trend and is essentially equivalent to 
the minimum mass for a SN II-P progenitor (see above; 
ISmartt et all [20091. It is notable that the case of SN 
2008bk offers a possible indication that me tallicity could 
also p lay a role in the low luminosity. ISmartt et all 
( 2009) assume in the cases of SNe 1999br and 2003gd 
that the metallicity in these environments was subsolar; 
on the other hand, they adopt essentially solar metallic- 
ity in the case of SN 2005cs. It will be interesting to 
see in future cases whether a two-parameter (mass and 
metallicity) trend is possibly emerging. 

We have already noted that SN 2008bk is ~ 1.0 mag 
more luminous than SN 1999br (My = —13.8 mag; 
iPastorello et all 120041) . This is particularly surprising, 
given the close overall agreement between the spectra, 
light-curve shapes, and color evolution for both SNe 
1999br and 2008bk; one is tempted to expect their lumi- 
nosities near maximum brightness to at least be compa- 
rable, if not the same (however, see lPastorello et alll2004l 
for an assessment of the diversity among low-luminosity 
SNe II-P). One possible reconciliation of the different 
luminosities could be to assume a distance to the host 
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galaxy (NGC 4900) of SN 1999br w hich is different 
from w hat has been adopted previously. iPastorello et al.l 
(2004) assumed a "short" distance (17.3 Mpc), which is 
essentially t hat derived from the Tully-Fisher relation 
(jTullvl 11988(1 : i f one adopt s the "long" distance (39.5 
Mpc) found by Uones et al.l (|2009f l using the Expanding 
Photosphere Method (EPM), then = -15.6 mag, 
which is now ~ 0.8 mag more luminous than SN 2008bk. 
Although EPM has its own inherent uncertainties and 
pitfalls, the implication is that, assuming SNe 2008bk 
and 1999br are very similar events, the actual distance 
to the SN 1999br host may be somewhere in between 
the long and short estimates. (Our by-eye adjustment of 
the SN 1999br light curves in Figure [5] to match those of 
SN 2008bk would then imply that the distance to NGC 
4900 could be ~ 32 Mpc.) 

Finally, although it will be ultimately satisfying to val- 
idate the identification of the progenitor star by obtain- 
ing images of SN 2008bk years hence, when it ha s faded 
close to invisibility (e.g., iMaund fc Smarttl |2009?Pl the 
quality of the pre-SN imaging data and high astromet- 
ric precision of the SN and progenitor positions instill in 
us confidence in the star's identification. Furthermore, 
the late-time / image from 2009 August (see Tabic [2]), 
shown in Figure I12[ when compared directly with Fig- 
ure [SJ demonstrates wit h relatively litt l e dou bt that the 
progenitor is known. As IMattila et a l. (2008) point out, 
in the i^-band pre-SN image, once the PSF of the RSG 
was subtracted away, a faint source ~ 0'.'5 south of the SN 
position remains (we also see this same source, when the 
PSF of the SN is subtracted from the late-time J-band 
images). Based on the arguments above, we consider it 
highly unlikely that this much fainter object is the SN 
progenitor. 
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TABLE 1 

Photometric Sequence around SN 2008BK a 



Star 


B (mag) 


V (mag) 


R (mag) 


/ (mag) 


A 


17.459(004) 


16.328(003) 


15.728(004) 


15.085(008) 


B 


16.939(003) 


16.283(003) 


15.910(003) 


15.447(008) 


C 


18.894(012) 


18.120(015) 


17.662(015) 


17.132(013) 


D 


18.574(005) 


17.898(003) 


17.512(004) 


17.045(008) 


E 


18.184(005) 


17.314(003) 


16.813(003) 


16.290(010) 


F 


19.415(007) 


18.251(003) 


17.531(003) 


16.832(008) 


G 


19.025(006) 


17.714(003) 


16.836(003) 


15.861(008) 


H 


20.092(011) 


18.642(004) 


17.664(000) 


16.561(008) 


I 


19.392(008) 


18.054(003) 


17.178(003) 


16.219(006) 


.1 


17.512(012) 


17.401(012) 


17.289(012) 


16.918(011) 


K 


15.013(001) 


14.278(018) 


15.080(051) 


13.497(014) 


L 


16.794(004) 


16.130(004) 


15.790(004) 


15.364(007) 


M 


17.515(004) 


17.060(004) 


16.814(004) 


16.444(007) 


N 


18.730(004) 


18.031(004) 


17.649(004) 


17.165(004) 



a Uncertainties (lcr) are given in parentheses, in units of 0.001 mag. 



TABLE 2 

Optical Photometry of SN 2008BK a 



UT date 


JD-2450000 


B (mag) 


V (mag) 


R (mag) 


/ (mag) 


Sourcc b 


2008 Mar 28.41 


4553.90 


12.94(13) 








CHASE 


2008 Mar 29.40 


4554.90 


12.89(14) 


12.99(10) 






CHASE 


2008 Apr 13.41 


4569.91 


12.94(05) 






LCO duPont 


2008 May 10.79 


4596.29 




12.99(13) 






S. Howerton 


2008 May 22.80 


4608.30 




13.02(06) 






S. Howerton 


2008 May 30.85 


4616.35 




12.97(26) 






S. Howerton 


2008 Jul 4.80 


4651.30 




13.07(12) 






S. Howerton 


2008 Jul 21.65 


4668.15 




13.52(13) 






S. Howerton 


2008 Aug 13.63 


4691.13 




15.64(11) 






S. Howerton 


2008 Sep 13.02 


4723.52 


17.70(07) 


16.02(06) 


15.11(05) 


14.46(08) 


Calar Alto 


2008 Sep 30.93 


4740.43 


17.83(07) 


16.22(06) 


15.32(05) 


14.62(08) 


Calar Alto 


2008 Oct 15.91 


4755.41 


18.08(08) 


16.43(06) 


15.56(05) 


14.84(08) 


Calar Alto 


2008 Dec 5.84 


4806.34 


18.38(07) 


16.87(06) 


16.06(05) 


15.40(08) 


NOT 


2009 Jul 02.37 


5014.87 




19.04(03) 


18.39(03) 


18.07(03) 


SMARTS 


2009 Jul 16.36 


5028.86 




19.21(05) 


18.52(07) 


18.22(03) 


SMARTS 


2009 Aug 13.13 


5056.63 


20.67(06) 


19.55(08) 


19.14(05) 


18.67(08) 


Calar Alto 


2009 Aug 13.16 


5056.66 




19.65(05) 


19.04(03) 


18.64(03) 


NOT 



a Uncertainties (lcr) are given in parentheses, in units of 0.01 mag. 

b CHASE = 0.41 m PROMPT2 + Apogee U47; LCO duPont = Las Campanas Observatory 2.5 m duPont + SITe2K; S. Howerton 
= GRAS-008 + SBIG ST-L-11K 3, GRAS-015 + SBIG ST-8 3 or SBIG ST-10 3 CCD; Calar Alto = 2.2 m Calar Alto + CAFOS; 
NOT = 2.5 m Nordic Optical Telescope + ALFOSC; SMARTS = 1.3 m SMARTS + ANDICAM. 



TABLE 3 

Optical Spectroscopy of SN 2008bk 



UT Date JD-2450000 Day a Grism/Grating Wavelength Instrument 13 

Range (A) 



2008 Jul 23.36 


4670.86 


20 


KPGL3-1 


4005-7738 


CTIO 


2008 Jul 24.35 


4671.85 


21 


KPGL3-1 


4005-7738 


CTIO 


2008 Jul 25.31 


4672.81 


22 


KPGL3-1 


4005-7738 


CTIO 


2008 Sep 14.04 


4723.54 


72 


G200 


4900-10458 


CAHA 


2008 Sep 16.02 


4725.52 


74 


G200 


3665-10463 


CAHA 


2008 Sep 30.93 


4740.43 


89 


R100 


5710-9585 


CAHA 



a The approximate day is referenced to our estimate of the end of the plateau phase, 
JD » 2454651.3. 

b CTIO = Blanco 4 m telescope + R-C Spec (Cerro Tololo, Chile); CAHA = Calar 
Alto 2.2 m telescope + CAFOS SiTc (near Almeria, Spain). 
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TABLE 4 

Photometry of the 
SN 2008BK 
Progenitor 51 



Band Magnitude 



V 22.81(09) 

R 22.03(07) 

I 20.71(08) 

J 19.26(14) 

H 18.55(06) 

K a 18.14(10) 



a Uncertainties (lc) 
are given in paren- 
theses, in units of 
0.01 mag. The JHK 3 
magnitudes are in the 
2MASS photometric 
system. 



TABLE 5 

Strong-line Indices and Oxygen Abundances 3 - for Hii Regions Near SN 2008bk 



Region SN offset (") b SN offset (pc) c N 2 12 + log(0/H) 03iV2 12 + log(0/H) Source 

A 12.3 203 -0.71 8.45 This work 

B 95.6 1576 -0.62 8.52 ■ ■ • • • ■ This work 

#27 82.2 1355 -0.66 8.53 0.81 8.47 Bibbv fc Crowther (20101 

W13 103.2 1701 <-0.77 <8.41 >1.10 <8.38 McCall, Rybski, & Shields (1985) 



a Computed following IPettini & Page] 120041) . 

b Angular dist ance of the H II region fro m SN 2008bk, deprojected for the host galaxy's position angle and inclination (adopting the values 

for these from lBibbv & Growtherll2010lK See Figure 

c Assuming a distance to the host galaxy of 3.4 Mpc (Pictrzvnski ct al. 2010). 




Fig. 1.— Chart of SN 2008bk in NGC 7793 in an /-band image obtained with the Calar Alto 2.2 m telescope on 2008 September 13. 
The comparison stars from Table[T] as well as the SN, are labeled. North is up and east is to the left. 
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i i I i i i I i i i I 

200 400 

Day Since Plateau End 

Fig. 2. — The limited BVRI light curves of SN 2008bk from photo metry in Table l2l (solid points). Shown for comparison are the light 
curves for the subluminous SN II-P 1999 br in NGC 4900 (open p oints; Pastorcllo ct al. 2004), a s well as the V light curves for the SNe II-P 
1999em in NGC 1637 (dot-dashed line; IHamuv et all 120011 : ILeonard et al.M2002t lElmhamdi etHTl I2003T ) and 2005cs in M51 (dashed line; 
Pastorcllo ct al. 2009). These latter comparison light curves have all been adjusted, by eye, in brightness to match the end of the plateau 
phase for SN 2008bk. All light curves are shown relative to the epoch corresponding to the end of the plateau phase. 
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Fig. 3. — The evolution of the absolute V magnitude of SN 2008b k (solid points). Sh own for comparison are the absolute V light curves 
for the sublumin ous SN II-P 1999br in NGC 4900 (open points ; Pastorello et al. 2003), as well as for the SNe II- P 1999em in NGC 163 7 
(dot-dashed line: IHamuv et al .1120011 ; I Leonard et al.1 12002; Elmhamdi et al. 2003) and 2005cs in M51 (dashed line: [Pastorello et al.ll2009T l. 
The observed V light curves for the SNe have been corrected for their respective extinction and distance; see text. All light curves are 
shown relative to the epoch corresponding to the end of the plateau phase. 
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Fig. 4.— Optical spectra of SN 2008bk obtained with the 4 m CTIO Blanco telescope on 2008 July 24, and with the 2.2 m Calar Alto 
telescope on 2008 September 16 and 30. See Table [3] and the text. Several notable features in the spectra are indicated. 
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4000 5000 6000 7000 

Rest Wavelength (A) 

Fig . 5. — A comparison of the SN 2008bk optical spectru m from 2008 July 24 (see Tabled with the spectra of SN 1999em l|Leonard et al.l 
2002) and a spectrum of SN 1999br (Pastorcllo ct al. 2004). The ages, in days, for the SNe in the figure are relative to the estimated epoch 
for the end of the photometric plateau phase in each case. 
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SN 2008bk 
, Progenitor 




Fig. 6. — The coadded i'-band image, from data obtained in 2007 September and October with the Gemini-South GMOS imager, showing 
the confirmed RSG progenitor of SN 2008bk (indicated by the arrow). The area of the image displayed here is 33" X 33". The effective 
seeing in the combined image is ~ 0'.'65. 
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HII regions near 
SN 2008bk 



W13 





Fig. 7 , — The continuum-su btracted Ho (+[Nll]) image of NGC 7793, an archival ancillary data product released by the Spitzer SINGS 
project (Kcnnicutt ct al. 2003), with the position of SN 2008bk indicated by a cross. The locations of Hn regions near the SN are labelled. 
Regions "A" and "B" wer e also detected in a s pectru m of SN 2008bk we o btained with the 2.2 m Cala r Alto telescope on 2008 September 
30. Region "#27" is from lBibbv fc Crowthe'rl ( feoToTl . Region W13 is from lMcCall. Rvbski. fc Shields! tl985ft . See Table |5] Note that the 
SN site appears to be within a faint, but extended, bubble of ionized gas, of diameter ~ 530 pc. North is up and east is to the left. 
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Fig. 8. — The observed spectral energy distribution (SED) of the SN 2008bk progenitor star {filled circles). For comparison, we also sh ow 
SEDs based on IRAF/SYNPHOT synthetic photometry extracted from the MARCS model stellar atmospheres (Gustafsson et al. 2008) at 
T c ff = 3600 K of 5 and 15 Mq stars (logg = 0.0) at solar metallicity ([M/H] = 0) with microturbulence velocities of 1—5 km s — 1 (dashed 
lines), and at subsolar metallicity ([M/H] = —0.25) w ith m icroturbulence velocities of 2 and 5 km s — 1 (dotted lines). We have reddened the 
model SEDs by our assumed value toward the SN (§ 14.11 1. The models were also normalized to the observed brightness of the progenitor 
in V. 
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Fig. 9. — An HRD showing the locus of the SN 2008bk progenitor. For comparison, we indicate model stellar evolutionary tracks, in 
panel (a) at metallicity Z = 0.02, in (b) at metallicity Z = 0.01, and in (c) at metallicity Z = 0.008. In (a), the 12 Mq models are from 
IHirsch i, Mcynct, & Macdcr (200 j|; "Geneva" models), with rotation ("rot," u; n j = 300 km s — 1 ; solid line) and without rotation ("norot;" 
long-dashed line), and from Eldridgc & Tout (2004;; "STARS" models), at initial masses 10 Mq (dot-dashed line) and 12 Mq (short-dashed 
line). In (b), we show just the STARS models at initial masses 7.5, 8, 8.5, and 9 Mq (dotted, long-dashed, solid, and short-dashed lines, 
respectively). In (c), we again show just the STARS models at initial masses 7.5, 8, and 8.5 Mq (dotted, long-dashed, and solid lines, 
respectively). 
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Fig. 9. — (Continued.) 
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Fig. 9. — (Continued.) 
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Fig. 10.— The observed SED of the SN 2008bk progenitor star, as measured by Mattila ct al. ( 2008; filled squares and upper limits). We 
show the SED we have measured, from Figu re [8] ( open circles), for comparison. We also display the templat e SED for a M4I supergiant in 
the LMC {Elias, Froecl, & Humphreys 1985; dotted line), extinguished by Ay = 1 ma g (Mattila ct al. 2008). Additionally, we illustrate a 
synthetic SED derived from MARCS model RSG stellar atmospheres (Gustafsson ct al. 2008) with T c g = 3500 K at subsolar metallicity 
([ M/H ] = —0.25), and with microturbulence velocities of 2 and 5 km s — 1 (dashed lines), extinguished by our assumed value toward the SN 
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Fig. 11. — An HRD showing the locus of the SN 2008bk progenitor, based on the A"-b and measurement by Mattila et al. (2008), and 
assuming the overall SED can be modeled by a MARCS model RSG stellar atmosphere (Gustafsson et al. 200 8) w ith T c g = 3500 K at 
subsolar metallicity ([M/H] = —0.25), and further ass uming our value of ext inction toward the SN (see Figure rro) . For comparison, we 
also sh ow the STARS massive-star evolutionary tracks (Eldridsc & Tout 2004), assuming a subsolar metallicity Z = 0.008 as Matti la" et al.l 
(2008) have, at initial masses 7 and 8 Mg (dashed line and solid line, respectively). 
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Fig. 12. — The /-band image, obtained on 2009 August 13 at the NOT, showing SN 2008bk (indicated by the arrow) at age ~ 507 d 
with / = 18.77 mag (see Tabled, only ~ 2 mag brighter than the progenitor shown in Figure [S] The scale and orientation here are 
approximately the same as those in the latter figure. The seeing was ~ l'.'O. 



